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ABSTRACT 

The  influence  of  surface  tension  on  solitary  and  cnoidal  waves  is 
discussed.  This  problem  was  first  considered  by  Korteweg  and  de  Vries  [1]. 
Shinbrot  [2]  attempted  to  derive  Korteweg  and  de  Vries  results  by  a  formal 
perturbation  expansion.  However  a  part  of  his  results  is  incorrect.  In  this 
note  an  approach  different  from  those  in  [1]  and  [2]  is  used  to  construct  a 
correct  perturbation  solution. 
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material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
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SIGNIFICANCE  AND  EXPLANATION 

In  this  note  we  present  a  perturbation  solution  for  gravity  capillary 
waves  in  water  of  finite  depth.  The  calculation  is  valid  when  both  the 
amplitude  of  the  waves  and  the  depth  of  the  water  are  small.  This  problem  was 
first  considered  by  KOrteweg  and  de  Vries  [1].  These  authors  showed  that 
waves  of  small  amplitude  in  relatively  shallow  water  can  be  described 
approximately  by  a  nonlinear  differential  equation,  in  addition  Korteweg  and 
de  Vries  obtained  periodic  solutions  of  this  equation  in  closed  form.  They 
named  these  solutions  cnoidal  waves.  Solitary  waves  are  the  limit  of  cnoidal 
waves  as  the  wavelength  tends  to  infinity. 

Shinbrot  [2]  attempted  to  derive  KOrteweg  and  de  Vries  results  by  a 
formal  perturbation  expansion.  However,  a  part  of  his  results  is  incorrect, 
because  the  possibility  of  a  depression  solitary  wave  was  excluded  by  him. 

In  this  note  an  approach  different  from  those  in  [1]  and  [2]  is  used  to 
construct  a  perturbation  solution  which  allows  depression  waves.  The  results 
obtained  agree  with  those  presented  by  Korteweg  and  de  Vries  [ 2] . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  MOTE  ON  SOLITARY  AND  CNOIDAL  WAVES  WITH  SURFACE  TENSION 
Jean-Marc  Vanden-Broeck  and  M.  C.  Shan 

§ 1 •  INTRODUCTION 

In  the  present  work  we  discuss  the  Influence  of  surface  tension  on 
solitary  and  cnoidal  waves.  Solitary  waves  with  surface  tension  were  first 
considered  by  Korteweg  and  de  Vries  [1].  They  showed  that  depression  solitary 
waves  exist  for  sufficiently  large  values  of  the  surface  tension.  Shinbrot 
[2]  attempted  to  study  solitary  waves  with  surface  tension  by  a  formal 
perturbation  expansion.  However,  a  part  of  his  results  is  incorrect,  because 
the  possibility  of  a  depression  wave  was  excluded  by  him.  In  this  note  we  use 
an  approach  different  from  those  in  [1],[2]  to  construct  a  perturbation 
solution  which  allows  depression  waves.  In  addition,  the  effect  of  surface 
tension  on  cnoidal  waves  is  also  considered.  The  results  obtained  agree  with 
those  presented  by  Korteweg  and  de  Vries  [1]. 

The  problem  is  formulated  in  Section  2.  The  perturbation  expansion  for 
solitary  waves  is  derived  in  Section  3.  Cnoidal  waves  with  surface  tension 
are  considered  in  Section  4. 

§2.  FORMULATION  OF  THE  PROBLEM 

The  governing  equations  for  the  steady  flow  of  an  incompressible  inviscid 
fluid  of  constant  density  with  surface  tension  in  reference  to  a  moving 
coordinate  system  at  a  constant  speed  c  in  the  x-direction  are  the 
following: 

u*  +  v*  -  0  ,  (1) 

*  y 

p(u*u*  +  v*u*)  *  -p*  ,  (2) 

x  y*x 
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X  y 

subject  to  the  boundary  conditions 


?t  +  u*Cx  "  v*  "  0 


>*  -  pg  , 

(3) 

y 

at  y  «  C*(x) 

(4) 

at  y  -  C*(x) 

(5) 

1  0  . 

(6) 

p*  "  “T*c£x(1  +  c$2>3/2  at  y  “  ?*{x)  (5) 

v*  -  0  at  y  -  0  .  (6) 

Here  (u*,v*)  is  the  velocity,  P  the  constant  density,  p  the  pressure, 

9  the  constant  gravitational  acceleration,  y  ■  C*  the  equation  of  the  free 
surface  and  T*  the  constant  surface  tension  coefficient.  Following  Peters 
and  Stoker  ([3]),  we  use  x  and  the  stream  function  t  as  independent 
variables,  and  transform  the  fluid  domain  to  a  fixed  horizontal  strip  in  the 
x,+-plane.  Furthermore,  it  is  also  more  convenient  to  use  u* ,  p*  and  the 
so-called  stream  line  function  f *  as  dependent  variables,  where  f* 
satisfies  4>(x,f*)  -  constant  to  define  a  stream  line.  If  we  introduce  the 
nondimensional  variables 

5  -  */h,  n  *  ♦/ch,  u  «  u*/c,  v  *  v*/c  , 

f  -  f*/h,  p  ■  p*/Pc2,  X  -  gh/c2,  T  ■  T*/Pgh2  . 

where  h  is  the  equilibrium  fluid  depth,  (1)  to  (6)  become 

u ?  -  -fnP5  +  f?Pn  ,  <  5  <  •  ,  (7) 

uf^^  +  u^f^  -  -Xf ^  -  Pn  ,  o  <  n  <  1  ,  (8) 

ttfn  -  1  ,  (9) 

f  »  0,  at  H  ■  0,  -•<?<•,  (10) 

p  -  -Tf^/d  +  f£)3/2  at  D  *  I,  -•<(<•.  (11) 

Assume  X  is  close  to  some  critical  value  A  and  let 

e  -  A  -  X,  o  -  5(e) 1/2  ,  (12) 

then  ( 7 )  to  (11)  in  terms  of  o  become 

ua  ‘  ‘Vo  *  Vi  '  (,3) 
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K 


e(ufoo  +  Vo1  “  efn  ■  Afn  “  pn 


uf  “  1  . 
n 


f  ■  o,  at  n  ■  o  , 


-Tef  /( i  +  eff)3/2  at  n  -  1  . 


§3.  SOLITARY  WAVE  SOLUTIONS 


To  solve  (13)  to  (17),  we  assume 


where 


u  "  u0  +  EU1  +  6  u2  + 


f  -  fn  +  ef „  +  e  f „  +  ... 
0  1  2 


P  ■  P0  +  eP1  +  e  p2  +  ... 


u_  »  1,  fft  ■  n,  pft  *  -*(o  -  i)  • 


Substitution  of  (18)  to  (21)  In  (13)  to  (17)  will  yield  a  sequence  of 


equations  and  boundary  conditions  for  the  successive  approximations.  The 
equations  for  the  first  approximations  are 

ula  ”  "  plo  “  *fi0  '  1 

1  ‘  “in  *  Pin  '  0  ’  1 

u,  .  t,0  -  °  .  ( 


From  (21)  to  (22),  we  have 


f ,(0.0)  -  0  , 


p,  (0,0)  -  0  . 


f  «  o  • 
iowi 


By  integration  and  making  use  of  (25),  we  obtain 


f,0  -  a’(o)n  , 

where  a,{0)  is  to  be  determined.  By  integration  again. 
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f1  *  a^ojn  +  b^n)  . 

Here  we  assume  that  a^(0)  ♦  0  as  0  +  •  or  a^(0)  is  a  periodic  function 

It 

with  period  L  and  /  a  (0)do  “  0.  It  follows  that 

0  1 

-  a^fljil  .  (28) 

From  (23),  (24)  and  (26),  we  have 


U1  "  “a1(0)  ' 


p1  -  [1  -  la ^0)1  (n  -  1)  . 

By  substituting  (28)  to  (30)  in  (22),  it  is  found  that 


If  a'(0)  f  0,  then 


a* (0) (t  -  1)  -  0  . 


1  . 


(29) 

(30) 

(31) 

(32) 


The  equations  for  the  second  approximations  are 

u20  “  “  p20  +  ai(0)n  "  •1(®)«5C®)  ,  (33) 

•|(0)n  -  a ,(0)  (34) 

f ^  +  u2  -  a2(0)  -  0  ,  (35) 

f2(0,O)  -  0  ,  (36) 

P2  -  -Ta^(0)  at  n  -  1  ,  (37) 


where  (28)  to  (30)  and  (32)  have  been  used.  As  before  we  find  that 


f 2  *  a2(0)n  -  a^(0)H3/6  ,  (38) 

u2  -  -a2(0)  +  a2(0)  +  a^(0)n2/2  ,  (39) 

P2  -  -a  (0)<n  -  1)  +  a  (o)(n  -  1)  +  a^(0)(-T  +  n3/6  -  n2/2  +  1/3)  .  (40) 

Substitution  of  (38)  and  (40)  in  (33)  yields 

(1/3  -  T)a™(0)  -  a'(0)  -  3*j(0Ua(0)  .  (41) 

If  we  impose  the  condition  a^(0),  a|J(0)  +  0  as  o  ♦  •  and  a^(0)  »  0,  then 
by  integration  we  obtain 
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rr 

>! 


9 


9 


i 


a,(0)  -  sech2(0/2)(1/3  -  T)“1/2 
It  follows  from  (19),  (21),  (28)  and  (42)  that 

f*  ~  h  +  h(1  -  A )sech2(x/2h) [ ( 1  -  X )/( 1/3  -  T)]1/2  . 

If  T  <  1/3  and  A  <  1,  there  is  a  solitary  wave  of  elevation.  In  this 
case,  c2  >  gh  and  the  solitary  wave  moves  at  a  supercritical  speed.  If 
T  >  1/3  and  A  >  1,  there  is  a  solitary  wave  of  depression  moving  at  a 
subcritical  speed  units  c2  <  gh.  We  may  define 

a  -  |h( 1  -  A )  |  , 

and  (43)  in  terms  of  a  becomes 

2  1/2 
f*  ~  h  ±  a  sech  ix/2h) [±a/( 1/3  -  T)]  '  . 

If  we  consider  a  periodic  solution  of  (41),  then  by  integration 

(1/3  -  T)[a'(0)]2  -  -a2(o)  +  aj(o)  +  c,a,(0)  +  c2  , 


(42) 


(43) 


(44) 


(45) 


where  c,,  c ^  are  two  arbitrary  constants.  The  periodic  solutions  are  called 
cnoidal  waves  and  will  be  discussed  in  the  next  section.  Finally  if 
T  -  1/3,  then  (41)  becomes 


a'(0)(1  -  38,(0))  -  0  , 


(46) 


and  the  only  possible  solution  is  a,(0)  ■  0. 


§4.  CNOIDAL  WAVE  SOLUTION 

For  a  bounded,  periodic  solution  of  (45),  it  is  necessary  that  the  right 
hand  side  of  (45)  must  have  three  real  simple  zeros,  and  the  solution  will 
occillate  periodically  between  two  of  the  zeros.  Therefore,  we  rewrite  (45) 
in  the  form 

(1/3  -  T)  [aj  (0 )  ]  2  -  [h,  -  8,(0))  [a,  (o)  -  h^  [a,(0)  -  Y)  ,  (47) 

where  h,,  h ..  are  the  extremum  values  of  a,(o)  and 


Y  -  1  -  h2  -  h, 


(48) 
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J 


Without  loss  of  generality  we  shall  assume  that  hj  and  hj  satisfy  the 
inequality 

(j  -  t)(1  -  X)(h1  -  h2)  >  0 

Following  Korteweg  and  de  Vries  [1]  we  introduce  the  parameter  X  and 

write 

a^o)  ■  h^cos  2X  +  h^sin^X  .  (49) 

Substituting  (49)  into  (47),  we  obtain 

6  ~  -  (1  -  K^sin^X) ^2  .  (50) 

aO 

Here  3  and  K  are  defined  by 


„  _  (A  1  -  3T-*1/2 

3  lT  h^yJ 
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(51) 


(52) 


We  take  the  origin  of  a  at  a  point  where  a^O)  attains  the  value 
h^«  Thus  integrating  (50)  we  obtain 


a 


(1-KCln^)1/2 


0F(X,K)  . 


(53) 


We  can  rewrite  (53)  in  terms  of  Jacobian  elliptic  functions  as 

cos  X  ■  cn  ^  .  (54) 

Substituting  (56)  into  (49)  we  have 

a,(o)  -  h„  +  (h.  -  h  Jen2  §  .  (55) 

1  2  1  4  P 

It  follows  from  (19),  (21),  (28)  and  (55)  that 


f*  ~  h  +  (1-X)hh2  +  (l-Xjhth^  h 2) cn2(^) [ ( h  1—  Y)(1-X)/(1/3  -  T)]1/2 


(56) 


The  wavelength  L  of  the  wave  Is  given  by 


1  •  ^30.,  i;  ”  -x>i  VSr,(,t>  -  i57> 

where 

TAX)  -  /  (1  -  K2sin^t)"1/2dX  .  (58) 

0 


The  detemination  of  the  solutin  is  completed  by  imposing  the  condition 
L 

/  8.(0) da  -  o  where  the  integral  is  over  a  wavelength*  Using  (49),  (50) 

0 

and  (52)  we  have 


0 


n/2  h^os^X  +  h^sin^X 
0  (1  -  X2sin^C)^2 


-  S&tYF^K)  +  (h1  -  Y  }E^(K)]  .  (59) 

Here 

B  (X)  «/  (1  -  X2sin^C)1/2dX  . 

0 

Hence 

YF^X)  +  (h1  -  Y)B1(K)  -  0  .  (60) 

The  eight  quantities  X,  h^,  h2,  Y,  X,  L,  T  and  B  are  related  by  the 
five  equations  (48),  (51),  (52),  (57)  and  (60).  If  for  instance  X,  L  and 
T  are  given  the  others  can  be  found. 

If  L  «•,  relations  (57),  (52)  (60)  and  (48)  yield  X  -  1,  Y  -  h2  -  0 
and  hj  ■  1.  The  cnoidal  wave  reduces  then  to  the  solitary  wave  considered  in 
Section  3. 

For  infinitesimal  amplitude,  equation  (41)  reduces  to 

(1/3  -  T)a*(o)  -  aj (0)  .  (61) 

The  solution  is 


(62) 


(63) 

(64) 

(65) 

(66) 


c2  -  ^  tanh  ieht1  +  —  k2)  •  (67) 

k  P  g 

In  the  limit  as  <h  ♦  0,  (67)  may  be  expanded  as 

1  +  (T  -  1/3)K^i2  +  odc^j2)  .  (68) 

gh 

A  comparison  of  (67)  and  (68)  indicates  that  the  cnoidal  wave  solution 
overlaps  the  classical  linear  theory  in  the  limit  as  the  amplitude  of  the 
waves  tends  to  zero. 
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